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Abstract 

We calculate the relic density of very heavy, stable scalar neutrinos in the minimal super- 
symmetric štandard model. We include all two-body final states, as well as the effects of 
co-annihilation with scalar electrons. We find that the sneutrino relic density is in the cos- 
mologically interesting región 0.1 < íl~h 2 < 1.0 for 550 GeV < m~ < 2300 GeV. For nominal 
values of the parameters of our galactic halo, recent results from the Heideľberg-Moscow 
direct detection experiment set an upper limit on íl~ which is lower by a factor ranging from 
two to ten, depending on m~. 



A sneutrino which is the lightest supersymmetric particle (LSP) is a theoretically attrac- 
tive dark matter candidate 0. Previous studies of sneutrino relic densities ^ have 
been concerned with a sneutrino with a mass m~ that is less than mw- Heavier sneutrinos 
can annihilate into a much greater set of final states; for m~ ~ 100 GeV, the annihilation 
rate is very large, and the relic density of sneutrinos is much too small to be of cosmological 
or astrophysical interest. However, as the sneutrino mass is increased, the annihilation rate 
drops and the relic density rises. We show that heavy (m~ > 500 GeV) sneutrinos may have a 
cosmologically interesting relic density. For nominal values of the parameters of our galactic 
halo, direct detection experiments [§] _ o already plače constraints on this range of sneutrino 
masses, and recent results from the Heidelberg-Moscow direct detection experiment || may 
even rule out the values of Q~ = p^/p C ňt that we find. 

The calculation of relic densities of heavy sneutrinos involves an additional complication 
not present for light sneutrinos, arising from the need to consider coannihilations || with 



selectrons. The masses of the sneutrino v and left-handed selectron e T , are related via 10 



m~ = m~ — m| cos 2(3 cos 2 6*w , (1' 



where tan/? is the ratio of the vacuum expectation values of the two Higgs doublets. We 
see that for tan/3 > 1 we have m~ L > m~, as we require if the sneutrino is to be the LSP. 
However, for m~ ^> mw, the mass splitting between v and e L will be small. Consequently 
their number densities are approximately equal until well after the sneutrinos go out of 
chemical equilbrium with all other particles, and we must follow the selectron density as 
accurately as the sneutrino density. 

We compute the tree level cross-section for annihilation to all allowed two particle final 
states. Table 1 lists the set of final states we consider. The f d and f u refer to T 3 = —1/2 and 
T 3 = +1/2 fermions, respectively. As intermediate states we allow sneutrinos, selectrons, 
the gauge bosons 7, and Z°, the three Higgs bosons of the minimal supersymmetric 
štandard model (MSSM), H? (i = 1,2,3), the char ginos Xi i} = 1;2), and the neutralinos 
X° (i = 1,2,3,4). There are also four-scalar point interactions. The couplings are those of 
the MSSM | rn)| , |TT| |. For simplicity, we ignore sfermion (in particular selectron) mixing, and 
for convenience we use the relation M\ = |M 2 tan 2 ^ w between the U(l) gaugino mass M x 
and the S U (2) gaugino mass M 2 which follows from grand unification. At tree level, the 
masses of the scalar Higgs particles H® and H® can be written in terms of mz, tan (3, and the 
mass the pseudoscalar H®. (Although the one-loop corrections to the Higgs masses can be 
large, they will not impact our results significantly, and so we do not include them.) We are 
then left with five free parameters: tan/3, m H o, m~, M2, and the Higgs mixing parameter /i. 
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Since our results are largely insensitive to tan /3, we fix tan (3 = 2 in our numerical work. 

Table 1: Initial and final states 



Initial State 


Final States 


vv* 


//, zz, w+w-, 

ZH^^W+H+^+H-, 
fl?flJ ) (i = l,2 ) 3) ) Ä?fl§, H^ 2) H° 3 


e L e* 


same as above, plus 77 , 7 Z 


ve* L 

vv 

e L e L 
ve h 


vv 

ee 
ve 



To derive a thermally averaged cross section, we make use of tlie technique of ref. [pL2 1 . 
We expand (crv re i) in a Taylor expansion in powers of x = T/m~: 

(av Te i) = a + bx + 0(x 2 ) . (2) 



Repeating tlie analysis of ref. [O for initial particles with different masses m\ and m 2 yields 



/ \ 1 fi 3(mi + m 2 )T\ 2 

(o-Vrei) = 1 w(s) s ^ (mi+m2)2+3(mi+m2)T + 0(T ). (3) 

Here 

w(s) = ^ / ^LIPS |.M| 2 (4) 
n } 1 ácos^cM |-M| 2 , (5) 



32tt S V2 y_ 

where ďLIPS is the Lorentz Invariant Phase Space element, p(s) is the magnitude of the 
three-momentum of one of tlie initial particles in the center-of-mass frame as a function of 
tlie total center-of-mass energy squared s, 9qm is the center-of-mass scattering angle, and 
\Ai\ 2 is the absolute square of the reduced matrix element for the annihilation, summed 
over final spins. The a and b coefficients may be read off the right-hand side of Eq. |3| after 
expanding in powers of x. 
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The sneutrinos and selectrons remain in chemical equilibrium with each other through 
freeze-out. We consider the total density n = J2i n i, where the index runs over u, v*,e h , and 
e*. We write the rate equation for n: 



where H is the Hubble parameter, r j = n^ q /n eq , nf 1 is the equilibrium density of particle 
species i, and <7y is the total cross-section for the i th particle to annihilate with the j th 
particle. Long after freeze-out, the selectrons decay into sneutrinos, leaving a number density 
n of sneutrinos. We rewrite the derivative in eq. (||) in terms of the temperature T and 
numerically integrate the resulting equation down to T = 2.75 K to get the relic number 
(and hence mass) density of sneutrinos. 

We find that the dominánt contribution to er cff comes from annihilations of uu*, e L e*, and 
z/e* to gauge bosons, and vv, e L e L , and ve h to neutrinos and electrons via neutralino and 
chargino exchange. 

In fig. 1 we show the relic density fl~h 2 as a function of sneutrino mass; h is the Hub- 
ble parameter today in units of lOOkm/s/Mpc. The dashed line assumes that two of the 
sneutrino generations are much heavier than the third, and that the lightest neutralino is 
just 10% heavier than the sneutrino. (We set /i = M 1 , but since the higgsinos do not couple 
significantly to the sneutrinos, the value of \x does not affect our results. Also, the pseu- 
doscalar Higgs mass was taken to be m H o = 500 GeV, but changing this parameter also 
has a negligible effect.) Note that fl~h 2 passes through the cosmologically interesting región 
0.1 < VL~h 2 < 1.0 for 550 GeV < m~ < 2300 GeV. Increasing M\ raises the masses of the 
gauginos and eventually turns off the gaugino-mediated annihilations of vv, e L e L , and Ve L to 
neutrinos and electrons. The dot-dashed line in fig. 1 is for this case, with M\ ~ 10 m~. We 
can also consider the case where all three generations of sneutrinos have very similar masses, 
and so are all in chemical equilibrium with each other at the time of freeze-out. Although 
one might expect that this case simply results in a value of íl~h 2 which is three times larger, 
in fact it is a little more complicated because some annihilation channels (Vfij — > ViVj, 
Vi(í Lj — > Vie j, e Li e Lj — > titj) can now occur in nine different ways, whereas all the others 
(ž/jž/* — > f f, etc.) can occur in three different ways. In fig. 1, the solid line corresponds to 
three equal-mass generations of LSP sneutrinos, but with the neutralinos again as light as 
possible (Mi — n ~ 1.10 m~). When the gauginos are heavy, then fl~h 2 is very close to what 
it would be for one generation of LSP sneutrinos. 



dn 



3Hn - (a eň v re \)(n 2 - n 2 ). 



(6) 



Hore 




(7) 
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The choice of masses for the sfermions here may seem a bit ad hoc. In the simplest grand 
unified model where all the sfermions share a common scalar mass parameter m, radiative 
corrections drive the mass of the sneutrino above the mass of the right-handed selectron [|13j . 
To keep the sneutrino as the lightest sfermion, and to prevent the sneutrino mass from being 
driven above the mass the lightest neutralino, we must invoke new high-energy physics to 
alter the running of the sneutrino mass. The GUT relation between M\ and M2 may then 
be affected, but we have seen that the effect of changing the gaugino masses by a factor of 
a few is small. 

Non-detection of halo particles in direct detection experiments have put bounds on the 
mass and relic abundance of dark matter sneutrinos ||, ||, ^ ||, ^ || . The dotted line in fig. 1 
marks the most stringent reported upper bound || on the relic density íl~ of sneutrinos for 
a given sneutrino mass m~. We assume the total dark matter density ÍŽdm — 1, so that the 
fraction of the dark matter halo which is composed of sneutrinos is íl~. As in ||, we take 
the local halo density of dark matter to be 0.3 GeV/cm 3 ; then = (0.3 GeV/cm 3 )fž~. The 
cross-section for sneutrino-nucleus interactions is four times greater than the cross-section 
for Dirac neutrino-nucleus interactions [13, [T5| (not two times greater, as sometimes stated 
0, U). In fig. 1, we take the published data of Jľ| to determine the upper bound on fl~. We 
compare their cross section constraint to four times the Dirac neutrino cross section. We 
see that there are upper bounds of m~ < 600 GeV and fl~h 2 < 0.08 GeV (for h = 1). Most 
cosmologically interesting relic densities are excluded. This assumes that the dark matter 
which is not composed of sneutrinos interacts less strongly with direct detection experiments 
than do sneutrinos. Also in fig. 1, the lower dotted line is found from íl = a^,/Aa c , where 
0"b and a c refer to the curves from fig. 3 of [§]. From this experiment we see that our 
computed range of íl~ is ruled out by a factor ranging from roughly two to ten (for h = 1). 
If we consider sneutrinos as the sole component of the dark matter, then we should fix 

= 0.3 GeV/cm 3 independent of Q~h 2 ; in this case we get a lower bound on the sneutrino 
mass of m~ > 17TeV, which we find requires a value of fl~h 2 that is much bigger than one. 
However, we should keep in mind that the uncertainty in the local value of our halo density 
is at least a factor of two fL6| . Furthermore, the momentum transfer in the interaction of a 
heavy sneutrino with a detector nucleus is large; a calculation of the interaction cross-section 
must therefore include nuclear form factors, which are not well known. 

To summarize, we have calculated the relic density of sneutrinos in a model where a 
heavy (> 500 GeV) sneutrino is the LSP. We have found that sneutrino masses in the range 
550 GeV < m~ < 2300 GeV may give cosmologically interesting relic densities, but that this 
range is ruled out for nominal values of the parameters of our galactic halo by a factor 
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ranging from roughly two to ten. Uncertainties in the halo parameters and nuclear form 
factors still leave a marginal chance for the heavy sneutrino to be a component of the halo 
dark matter. 
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Figúre Caption 



Fig. 1) The relic density fl~h 2 as a function of m~. The dashed line corresponds to 
m H° = 500 GeV and m~o ~ 1.10 m~, with one generation of LSP sneutri- 
nos. The dot-dashed line is for m~o ~ 10 m~. The solid line corresponds to 
the same parameters as the dashed line, but with three generations of LSP 
sneutrinos. The Neuchatel-CIT and the Heidelberg-Moscow [|J direct de- 
tection experiments exclude values of Q~ (not £l~h 2 ) above the upper and 
lower dotted lines respectively. 
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